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Inhomogeneous charge distributions have important repercussions on electrostatic interactions in systems of 
charged particles but are often difficult to examine theoretically. We investigate how electrostatic interactions 
are influenced by patchy charge distributions exhibiting certain point group symmetries. We derive a general 
form of the electrostatic interaction energy of two permeable, arbitrarily charged spherical shells in the Debye- 
Hiickel approximation and apply it to the case of particles with icosahedral, octahedral, and tetrahedral 
inhomogeneous charge distributions. We analyze in detail how charge distribution symmetry modifies the 
interaction energy and find that local charge inhomogeneities reduce the repulsion of two overall equally 
charged particles, while sufficient orientational variation in the charge distribution can turn the minimum 
interaction energy into an attraction. Additionally we show that larger patches and thus lower symmetries 
and wave numbers result in bigger attraction given the same variation. 



I. INTRODUCTION 

Most commonly used theoretical models of virus cap- 
sids usually consider charge to be homogeneously dis- 
tributed across one or two spherical shells represent- 
ing the viral capsid shelP^. Reality is more complex, 
and the full 3D resolution of capsid structure shows sig- 
nificant variation in the distribution of charged amino 
acids across and along the capsid, carrying a si gn ature 
of the underlying symmetries of the structure^. Lo- 
cal variations of charge on the capsid can play an im- 
portant r ole in different scenarios, from RNA-capsid in- 
teractions^!, interactions of the capsid with polyvalent 
ions^, to capsid-capsid electrostatic interactions^ or 
nanoparticle-templated assembly of virus-like particles^. 

In addition, electrostatic interactions can be the sta- 
bilizing force in ionic colloidal crystals of binary mix- 
tures of oppositely charged particles which have been 
observed both in experiment and in simulatiorP^HH] j n 
such cases the assumed interaction potential is usually 
taken to be spherically symmetric, facilitating assembly 
but also severely limiting potential symmetries of assem- 
bled crystals^. It is thus of special interest to inves- 
tigate more general forms of the interaction potential 
between colloidal particles that are often anisotropic in 
shape and c ontai n inhomogeneously distributed interact- 
ing nioietiePMl. An important example are the patchy 
colloids as well as inverse patchy colloids with hetero- 
geneously distributed interacting patches that can form 
a richer variety of admissible structures"^ Charged 
patchy colloids with heterogeneously charged surfaces be- 
long to the same variety of complex colloidal particles 
except that in charged systems the interaction is usually 



much longer ranged than the standard bonding distances 
in ordinary patchy systems with limited valence^. 

Contrary to classical colloids, virus capsids, which are 
the main motivation for our work, offer a number of ad- 
ditional features that colloids do not possess. A few more 
prominent properties are the unique, spatially-defined 
chemistries on their surfaces as well as absolute monodis- 
persity in size and mass distributions^. Such monodis- 
perse viral colloids are highly symmetrical^ 2 ^, with pro- 
nounced symmetry-related variation in particle surface 
chemistry. These surface structural features make them 
a possible comp onent material for the production of bio- 
photonic crystals^25|, non-close-packed crystalline struc- 
tures (see Ref. [IT] and references therein) and virus-like 
nano particles that promise to invigorate the search for 
the perfect gene-therapy vector^. 

In vitro the interviral interaction potential is charac- 
terized by weak long-range electrostatic repulsion and 
sterically mediated close-range attraction^. Our focus 
in this work will reside on the electrostatic part, and we 
will derive an analytical (closed-form) expression for the 
interaction free energy between two spherical shells carry- 
ing arbitrary inhomogeneous surface charge distributions 
within the linearized Poisson-Boltzmann (Debye-Hiickel) 
approximation 2 ^. This result will enable us to analyze 
models of highly symmetric charge distributions mimick- 
ing the symmetries of viruses within the ap proach re- 
cently proposed by Lorman and RochapsH^!] j n order to 
describe the mass distributions in spherical viruses. 

The rest of the paper is structured as follows: In order 
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to introduce the setting, we first summarize in Sec. II A 
results obtained in previous works for the electrostatic 
self-energy of a single arbitrarily charged spherical shell 
in the Debye-Hiickel (DH) approximation. Building on 
this we derive an analytical expression for the electro- 
static interaction of two such shells in Sec. |IIB| Ap- 
pendices [5] and [B] contain the details of the derivation 
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and certain limiting cases. We first use the derived ex- 
pression to analyze a simple model of axially symmet- 
ric quadrupole surface charge distribution in Sec. |III| 
comparing some of the results with those obtained in 
the literature fo r a similar case, but with a different 
modePM In Sec. IV we then introduce more general sur- 
face charge densities with tetrahedral, octahedral, and 
icosahedral symmetry and use them to examine in more 
detail the electrostatic interactions between shells carry- 
ing such charge distributions. We end with a discussion 
and conclusions in Sec. IVl 



II. DERIVATION OF INTERACTION ENERGY IN 
DEBYE-HUCKEL APPROXIMATION 

Throughout the paper we will deal with electrostat- 
ics of charged shells in a monovalent salt solution within 
the framewo rk of the linearized mean- field Debye-Hiickcl 
(DH) theorjl2Z12oH32]_ This approximation is reasonable 
for sufficiently small surface charge densities on the sur- 
face of the particles, low ion valencies, high medium di- 
electric constant, or high temperatures. All the require- 
ments are in general well fulfilled for the case of mono- 
valent salt solutions and the surface charge densities rel- 
evant for viruses or colloidal particles considered here. 



A. Self-energy of a charged spherical shell 

We first briefly sum up the relevant parts in the 
derivation of the electrostatic self-energy of an arbitrarily 
charged spherical shell in a salt solution in order to intro- 
duce some concepts important for the rest of the paper. 
The radius of the shell is R and the surface charge den- 
sity on the shell can be expanded in terms of spherical 
harmonics, 



a(Cl) =J2a(lm)Y lm (n). 



(1) 



The requirement that the surface charge density be real 
imposes the constraint a*(lm) = (— I)"V(^to), where we 
have introduced the notation m = —m. The electrostatic 
potential ip is given by the solution of the DH equation 



K 2 (p, 



(2) 



where k is the inverse DH screening length of a monova- 
lent I : 1 salt with bulk concentration cq, k = \/8ttIbCq, 
with Ib = eQ/AireeoksT the Bjerrum length and T the 
absolute temperature. The dielectric constant of water 
is taken to be e — 80. 

The electrostatic potentials inside and outside the 
shell, 



(pi(r, 0) = )^a(lm)ii(Kr)Yi m (Q), and (3) 
n) = J2 b (lm)k l (Kr)Y lm {n), (4) 



are expanded in terms of spherical harmonics Yi m (Cl) and 
the modified spherical Bessel functions of the first and 
second kincP^, 



il{x) 
h(x) 



4f l/aOO and 
Ki+i/ 2 {x). 



(5) 
(6) 



The potential has to be continuous on the surface of the 
shell, <pi(R, O) = (fn(R, CI) for all values of Cl, and the 
surface charge on the shell gives rise to a discontinuity of 
the electrostatic potential on its surface in the standard 
form 
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By applying the boundary conditions we find that the 
solutions for the potential have the same symmetry as 
the underlying surface charge density, 



<Pi(r,Q) 
tpu{r, Cl) 
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Here, we have defined 



keeq ki(nR) 



a(lm)Y lm (n), (8) 
a(lm)Y lm (ri). (9) 



C (l,x) = xI l+1/2 (x)K l+1/2 (x). 



(10) 



The electros tatic free energy in the DH limit can then be 
calculated as^U 



F DH = \i <j(Cl)p(R, Q)dS. 

1 JdV 



(11) 



Using the orthogonality relations for the spherical har- 
monics together with the requirement that the charge 
density be real, we then obtain 



F. 



R 2 



DH 



4^o 



2k££o 1 + coth nR 
;>o m 



\a(lm)\ 2 . (12) 



Lm 



This result was already obtained by Marzec and Day^; 
the first term (/ = 0) is simply the self-energy of a uni- 
formly charged shell in a salt solution^. Details of the 
derivation can be found in e.g. Ref. [3l where it was used 
to derive the free energy of a partially formed shell - a 
special case where the surface charge density is a Heavi- 
side step function in the azimuthal angle. 



1. Rotations of the shell(s) 

Given the original surface charge distribution of a shell 
[Eq. ([I])] in its reference frame, the expansion coefficients 
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will in general change with different orientations of the 
shell as we rotate it. The easiest way to incorporate the 
rotations is through the Wigner matrices fl2 ra (w) given 
in terms of the Euler angles u) = (a, /3, 7) in the zyz 
notation: 



(13) 



where we use the definition of the small d-matrices 
^m'ro(/3) as given by Rose^l. From here the expansion 
coefficients of the rotated dist ributi on can be expressed 
in terms of the original ones a d 37 * 38 ! 



(14) 



the rotations do not mix different wave numbers I. 
Due to the way the expansion coefficients feature in the 



self-energy of the shell [Eq. ( 12 )] changing the orientation 
of the shell does not influence the result. However, when 
we will next consider the interaction of two such shells, 
different orientations of both shells will play a role. In the 
rest of the paper we will thus denote the coefficients of 
the original surface charge distributions of the two shells 
in their respective reference frames by a(lm, i) and those 
of the rotated distributions by a'(lm,i), while the orien- 
tations of the two shells will be given by the Euler angles 



B. Interaction free energy of two arbitrarily charged 
spherical shells 

Taking two charged shells as studied in the previous 
Section, we now derive the expression for their interaction 
energy. The basic aspects of the method were already 
presented by Langbein^ in the calculation of van der 
Waals forces between molecules, and similar approaches 
have been used to calculate the force and interaction en- 
ergies o f charged shells / colloids in simpler^SHM anc j more 
involved 21 * 45 " ^ scenarios. 

The two shells under consideration carry surface charge 
densities a^\fl) expanded as in Eq. 0; the orientations 
of the two shells can be arbitrary and are given by the 
Euler angles u)j — (aj, The shells and their re- 

spective properties are denoted by indices 1 and 2. The 
radii of the shells are Rj = R, and the distance between 
their centers is r 2 i = p. The vector p = r 2 r x = (p, £, rj) 
in spherical coordinates connects the two coordinate sys- 
tems. The rest of the system properties are the same as 
in Sec. |II A| Fig. [T] shows a sketch of the system. 



1. Electrostatic potential . . . 

The solution of the DH equation for the mean elec- 
trostatic potential [Eq. ([2])] can be written as a linear 
superposition ip = ipi + ip2, where ipj is the electrostatic 
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Figure 1. The system under consideration: two shells with 
arbitrary surface charge densities a^\Cl) (in this case, with 
icosahedral and octahedral symmetry) and radii R are located 
at a distance p. Either or both of the distributions can be 
additionally rotated, with the orientations given by the Euler 
angles u>j. The shells are in a monovalent salt solution with 
bulk concentration Co giving rise to the electrostatic screening 
length k . The background color of the shells represents the 
homogeneous charge contribution, whereas the other tones 
depict the local patches of charge of opposing signs. 



potential of a single (j-th) shell. We can divide the sys- 
tem into three regions: inside the first shell {ipi n ,i)i inside 
the second shell (V>in,2)j and outside both shells ("J). The 
boundary conditions are again the continuity of the po- 
tential at the surface of each shell, 



i/j in , j (Rj,n) = *(R+,n), 



(15) 



and the discontinuity of the derivative, which has to be 
proportional to the charge on the shell, 



in,] 



dr-j 



(n) 



ee 



(16) 



The electrostatic potential can be expanded in terms of 
modified spherical Bessel functions of the first and sec- 
ond kind and spherical harmonics. We will follow the 
notation introduced in Ref. [48] and define 

i}j+ n {v) = n{Kr)Y lm {tt) and (17) 

^(r) = *,(«r)y Im (n). (18) 

Thus we can write for the solution inside each shell 



V'sraj ( r j) 



a{lm,j)i}'l m {v j ), 



(19) 



and the solution at a point r outside both shells 
*(r) : > '[6(Zm,l)^- m (r-n)+6(Zm,2)^" m (r-r 2 )] 



E 

l.m 



(20) 

due to the linearity of the DH equation. We will also use 
the notation ^^,(1) and tpf^i^) for the potentials written 
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in the coordinate systems of the first and the second shell, 
respectively. 

To move between the two coordinate systems we have 
to use an addition theorem for the solutions of the 
modified Hclmholtz equation in spherical coordinates^!. 
Clercx and Schram^ give the theorem for an arbitrary 
relative positioning of the two shells specified by the vec- 
tor p. However, since all possible orientations of the two 
shells are already taken into account via rotations given 
by the Euler angles uij we can fix their relative position 
to p = pi. This allows us to write a simplified version 
of the addition theorem from Ref. HH1 



(21) 



where, similarly as in the case of a single shell, 



Feij = l<f a^iSl^^r^R^dSi, (25) 

where now ip( r j — Ri ^j) designates the potential of both 
shells so that their electr ostati c interaction energy at a 
distance p is then given a d 43 * 49 ! 



V mt (p) = F el (p) - F el (oo). 



(26) 



First we note that the expansion coefficients a(lm, 2) 
[Eq. ( 24 )] can be written in the form 



where the function is expressed in terms of mod- 

ified Bessel functions of the second kind and Wigner 3-j 
symbols: 



cz (i) = ]T(-iy +m (2 S + iW(2i + i)(2p + i) x 

s 



The theorem in the opposite direction, ip lm (2) 
J2 P q (2)i/;+ q (l), is essentially the same, with an addi- 



tional factor of (—1) 
s in the function C^(2) 
of the function Cf^L are listed in Appendix 
Combining now 



appearing in the summation over 
Some of the useful properties 



Al 



the two boundary conditions in 
Eqs. (|15| and (|16| on the surface of e.g. the second shell, 



we obtain the coefficients of the solution for the potential 
in the exterior and interior, 



b(lm,2) 
a(lm, 2) 



a(lm,2) C q (1,kR) 

ksso ki(nR) 
a(lm,2) C (1,kR) 



E 

P,Q 



neso ii{nR) 
a(pq,l) C (p, kR) 
keeq kp(nR) 



(23) 



CW. ( 24 ) 



where the function Cq{1,x) is defined in Eq. (10 1. The 
expressions for the coefficients of the first shell are simi- 
lar with 1 4-» 2 exchanged everywhere. The procedure for 
obtaining the coefficients is given in Appendix | A 2| This 
completes the solution for the electrostatic potential in- 
side and outside the shells. 



2. ... and the interaction energy 

From the solution for the DH potential we can calculate 
the free energy of the two shells as F e i — F e n + F e i2 



a(lm, 2) = a (Zm, 2) + £ b(pq, 1)££(1), (27) 

with ao(lm,2) being the expansion coefficients one ob- 
tains in the case of a single shell (Sec. II A). By writing 
out in full both the surface charge distribution and the 
potential at the shell surface in Eq. ( 25 ) we can then use 



the properties of spherical harmonics to write the free 
energy of the second shell in the form 



p - p(°) 



R 2 



^^6( M ,l)a*(/m,2)4™(l) ii (^), 



(28) 

where Fj?\ is now the free energy in the case of a single, 



isolated shell [Eq. Q]. 

In the interaction energy, these terms cancel out with 
the contribution of the two shells at infinite separation, 
and we are left with 



V int {p) 



R 2 



££[&( OT ,lK(Zm,2)4™(l) 



l,m p,q 



+ b(pq,2)a*(lm,l)£ L ™(2) %i{kR) 



(29) 



We can also write this expression in terms separated by 
their respective wave numbers I and p, and define 



V in t = ^11 + ^2^+^0 = I 



2 £ 

I 



i>p 



(30) 

The functions W\ v can be further simplified using some of 
the properties of Wigner 3-j symbols and Bessel functions 
(see the Appendix A3), finally leading to the interaction 
energy terms 
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R 



W iv = C(l,p, kR) V(-l)' +m Re a(lm, l)a*(pm, 2) + {-l) l+p cr(pm, l)a*(lm, 2) 



x (2s + 1)^/(2/ + l)(2p+l)fe s (/cp) ( 



I s p\ I I s p 
I { to to 



(31) 



where we have defined 

C(l,p,x) = xl l+1/2 (x)l p+1/2 (x). (32) 

The entire dependence on the orientations of the two 
shells relative to their respective reference frames is hid- 
den in the expansion coefficients a(lm,i), and upon re- 
placing them with the rotated ones, a'(lm,i), we get the 
interaction free energy terms Wi p (p, U3\, u)2). 

III. SIMPLE CASE: AXIALLY SYMMETRIC 
QUADRUPOLE DISTRIBUTION 

With the expression for the interaction free energy 
known, we can now analyze different scenarios. Firstly, 
we will shortly discuss an example of two shells carrying 
axially symmetric quadrupole distributions, a case simi- 
lar to the one considered by Bianchi et al.^ for inverse 
patchy colloids. There, however, the colloids are imper- 
meable to salt ions, unlike the shells considered here. 

The model can be described by a surface charge distri- 
bution of the form 




p/2R p/2R 



Figure 2. Quadrupole-quadrupole interaction energy V22 as 
a function of the intershell separation for three different con- 
figurations of the shells, pole-pole (PP), equator-pole (EP), 
and equator-equator (EE). The radii of the shells are R = 10 
nm, with a) kR m 10 and b) kR sa 3, and the strength of the 
variation is <72 = 1 eo/nm 2 . 



a(ft) = <j 2 Y 2Q (Q) + V^a Y aQ (Q), 



(33) 



which is axially symmetric around the z axis and has 
two patches of the same charge at the poles and a patch 
of opposite charge around the equator. By varying <t 
we can also include a homogeneous background, where 
we set the coefficient <t(00, i) = \Tkno§ to ensure that 
the total charge is given by Q = 4irR 2 (To. Due to the 
symmetry the rotation angles are limited and we need 
consider only the azimuthal angles of the two shells, 

Focusing initially only on the quadrupole-quadrupole 
interaction energy V 22 (p, u>j) we can see from the symme- 
try of the problem that the most optimal arrangement of 
the shells in this case is equator-pole (EP). Thus, we 
consider three orientational configurations (EE, EP, PP) 
and plot the radial dependence of the interaction energy 
for the three cases in Fig. [2] Similar to observations in 
Ref. l2"T1we find that the interaction energy is most repul- 
sive in the PP configuration, is smaller for the EE config- 
uration, and turns into attraction in the EP case. With 
increasing kR the scale of the interaction energy dimin- 
ishes, and the interaction falls of more quickly with in- 
creasing inter-shell separation. The interaction between 
the shells also remains either purely repulsive or purely 
attractive with respect to the separation between them. 

By adding a homogeneous background to the two shells 
we obtain two additional contributions to the interaction 



energy, 

VintfaWi) = V {p) + W 2a {p,u}{) + V 22 {p,u,i). (34) 

The relative contributions of the terms are governed by 
the ratio <J 2 /ao, and Fig. [3] shows how the total inter- 
action energy changes with respect to this ratio for the 
EP configuration of the shells. The homogeneous con- 
tribution to surface charge density ctq = 0.1 eo/nm 2 is 
kept fixed and amounts to a charge of Q ~ 125 eo on a 
shell of radius R — 10 nm, while we vary the quadrupole 
component of charge on the shells. The homogeneous 
background (of the same sign on both shells) adds an 
additional repulsive contribution which prevails for small 
variations in surface charge density; thus, even the EP 
configuration of the shells results in a repulsive interac- 
tion. However, upon increasing the quadrupole part the 
interaction becomes attractive for larger and larger dis- 
tances, with the repulsive contribution slowly vanishing. 

Figure [4] shows the contribution of the different terms 
in the interaction free energy [Eq. (34)] for the PP, EE, 
and EP orientations of the two shells. The homoge- 
neous term Vq contributes only a repulsive component 
regardless of their orientations, whereas the cross-term 
W 2 q changes the behavior of the interaction energy with 
respect to the shell orientations. In the PP configuration 



() 




p/2R 

Figure 3. Total interaction energy of two charged shells with 
quadrupole variation [Eq. ( 34 1] as a function of the intershell 



distance for different ratios of the quadrupole and homoge- 
neous contributions to the surface charge density, o^/co 6 
[0.45,0.65]. Radii of the shells are R — 10 nm with kR ~ 3, 
and the shells are in the EP configuration. 



every component is repulsive, totalling to an even big- 
ger repulsion than in the simpler quadrupole-quadrupole 
case. In the EP configuration, the total interaction is 
attractive due to the quadrupole-quadrupole term, the 
cross-term being consistently repulsive. The cross-term 
becomes attractive in the EE configuration, where the 
quadrupole-quadrupole term is now repulsive; the lat- 
ter, together with the homogeneous contribution, prevails 
and the total interaction is repulsive. 



IV. INTERACTION BETWEEN TWO SHELLS WITH 
SYMMETRIC CHARGE DISTRIBUTIONS 

In this Section we will examine electrostatic interac- 
tions between model virus-like particles carrying surface 
charge distributions with icosahedral symmetry; we will 
also include octahedral and tetrahedral symmetries for 
comparison. In order to do this, we next show how to ob- 
tain distribution functions invariant under a given sym- 
metry group. 



A. Lower symmetry charge distributions: tetrahedral, 
octahedral, and icosahedral symmetry 

In constructing the surface charge distribution function 
possessing a given symmetry, we will use the ideas pre- 



sented in Refs. [25] and [521 with appropriately described 
details. The sought-for symmetry adapted functions are 
obtained by considering crystallization on a sphere with 
symmetry reduction from the (isotropic) 3D rotation 
group SO (3) to the icosahedral (octahedral, tetrahedral) 
point subgroup y C SO(3f®. 

The symmetry adapted icosahedral functions are thus 
obtained by reducing the irreducible representations D"> 
of SO (3) to the icosahedral point group y. We stipulate 
that the identical representation of the point group ap- 
pears at least once in the reduction so that the resulting 
structure is invariant under all symmetry operations of 
the icosahe dral g roup. This restricts the allowed wave 
numbers I tcP122J 



l lco = 6^ + lOj (+15). 



(35) 



Odd Is lack inversion symmetry, making such distri- 
butions consistent with the asymmetry of the proteins 
which constitute the viral capsicP^. Even Is are included 
for comparison and the corresponding distributions fit a 
dodecahedral geometry better than the icosahedral. For 
octahedral and tetrahedral groups the same procedure 
yields the allowed wave numbers l oct — Ai + 6j (+9) and 
l tet = U + 6j(+3). 

For a selected wave number the explicit form of the 
surface charge density <Jt(Sl) is then given by the basis 
functions f*(SY), k= 1, . . . , mx(l), of all itia totally sym- 
metric representations of the icosahedral group y in the 
restriction of the "active" irreducible representations of 
the SO (3). For lower wave numbers there is additionally 
only one identical representation for each allowed wave 
number, ttiaQ) = 1, and we have 



on(fi) = on/j(n). 



(36) 



The explicit form of the symmetry adapted irreducible 
icosahedral density function ft (SI) for a given wave num- 
ber is obtained by averaging the spherical harmonics over 
the icosahedral symmetry group. Thus we obtain dis- 
tributions of charge on a sphere which have icosahedral 
symmetry, characterized by an allowed wave number I: 



at(Sl) 



Q\ £ JlmYlm(Sl), 



(37) 



where the expansion coefficients are obtained with the 
averaging procedure as 



r = JL V D {1) 

l m en / < mm 

Gey 



,(G(a,/3, 7 )). (38) 



Here, G(a,j9, 7) are the symmetry operations of the 
icosahedral group given in terms of the Euler angles and 
Dmm' are tne Wigner matrices^. The irreducible den- 
sity functions for the octahedral and tetrahedral group 
are obtained in an analogous fashion. A few examples 
of the irreducible density functions /; of the symmetry 
groups are shown in Fig [5] 
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1 1.2 1.4 1.6 1 1.2 1.4 1.6 1 1.2 1.4 1.6 



p/2R p/2R p/2R 



Figure 4. Total interaction energy of two charged shells with quadrupole variation Vint (full line) and the relative contributions 
(Vb - dotted line, W20 — dot-dashed line, V22 - dashed line) as a function of the intershell separation for three different 
orientations. As before, the radii of the shells are R = 10 nm with kR « 3, with the ratio 02/00 = 0.55. 



It is important to stress that the expansion coefficients 
of the irreducible density functions /;(£!) are com- 
pletely determined by the required invariance under the 
point group and the only free parameter left is 07. The 
coefficients are purely real for even I and purely imagi- 
nary for odd with an additional useful property that 
S m \hm\ 2 = 1' i- e - tne functions fi are normalised. As 
can be seen in Fig. [5j the number of local minima and 
maxima increases with increasing wave number; for the 
icosahedral symmetry, their number can be correlated 
with the Caspar-Klug triangulation number (cf . Ref. |2"9"|) . 
Additionally, the functions with odd I have equal ex- 
trema, while those with even I differ in the depth and 
shape of the minima and maxima. 

In what follows we will deal with surface charge distri- 
butions of the form 



(39) 



where the wave number I and the symmetry fully deter- 
mine the expansion coefficients. The total charge on the 
shell is Q = 4irR 2 <ro and the 07 gives only the variation 
of the charge due to the symmetry, contributing nothing 



to the total charge (d 2 . 



of). We can 



connect the variation strength 07 with the symmetry, if 
we demand that the charge on the shell is of the same 
sign everywhere, 



a t fi(n) + (x > vo, 



(40) 



from where we obtain the critical variation strength 07 

_ / I fmin I 



B. Self-energy 



From Eq. ( 12 ) and the aforementioned property of the 



functions fi(fi) that J2m \^lm\ 2 — 1> we have for the self- 
energy of a single shell with surface charge distribution 
with wave number I 



7 (i) 

DH 



C {1,kR). 



(41) 



Even though the symmetry has an imprint on the func- 
tions /; and consequently on the electrostatic potential, 
the free energy depends only on the wave number and the 
strength of the variation. The function Cq(1,x) tends to 
1/2 in the limit when x — > 00, regardless of the wave 
number. However, in the opposite limit of x — >■ it tends 
to zero as x/(2l-\-l) in the lowest order of approximation. 

Upon the addition of homogeneously distributed back- 
ground charge, the self-energy of such a shell can be writ- 
ten as 



r (0) 

DH 



1 
47T 



<7,y c {i,kR) 

00) C (0,kR) 



(42) 



The ratio Cq{1,x)/Cq{Q, x) goes to 1 in the limit of x — > 
00, and falls off to an /-dependent constant in the limit 
of x — > as (l + x)/(2l + l). Consequently, in the limit of 
kR 3> 1 only the strength of the variation 07 and the ratio 
07 / Co play a role as the self-energy becomes insensitive to 
the wave number of the variation. In the limit of kR — > 
the self-energy with a given I goes to 0, the more so the 
larger the wave number, but with the scale still set by 



Figure 5. Symmetry adapted irreducible density functions fi(Q) for different symmetries and wave numbers: a) Icosahedral 
symmetry with I — 15, b) octahedral symmetry with I = 9, c) icosahedral symmetry with / = 21, and d) dodecahedral 
symmetry with / = 6. The distributions are mapped from a sphere to a plane using the Mollweide projection^. Warmer colors 
correspond to regions with positive values and colder colors correspond to regions with negative values. 



ai and (Jq. Figure [6] shows the free energy correction of 
a charged shell due to local charge variation as given by 
Eq. El, 



F DH = F < g> 1 (l + lcorr ), 



(43) 



in the limits discussed above and for the special choice 
of a; /do = l/\ff nin \. The local variation increases the 
shell self-energy, and quite significantly so in the limit 
kR ^> 1. On the other hand, the correction due to 
variation is smaller in the limit of kR <C 1, and for the 
range of wave numbers relevant for the icosahedral sym- 
metry is of the order of less than 1% of the homogeneous 
contribution. The correction decreases with increasing 
wave numbers as smaller "wavelengths" of the extrema 
become smoothed out due to screening. The scattering 
of the data with respect to the wave numbers is a con- 
sequence of the minima of the functions fj nm changing 
non-monotonically with I. 



C. Interaction of two equal shells 

Let us finally turn attention to the interaction of two 
equal shells carrying charge distributions with icosahe- 
dral, octahedral, and tetrahedral symmetries. We will 
be mainly concerned with determining the minima of the 
interaction free energy as a function of the orientation 
of the shells. Two cases of inhomogeneous charge dis- 
tributions will be explicitly considered: (i) the case of 
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Figure 6. Correction to the free energy of a charged shell 
[Eq. (43 1] due to local charge variation characterized by wave 
number The strength of the variation is such that the charge 
on the shell is everywhere of equal sign, oi — ao/\f[ rLm \. Full 
symbols show the limit where kR — > oo and empty symbols 
show the limit of kR — > 0. Circles denote icosahedral, squares 
octahedral, and diamonds tetrahedral symmetry. 



an overall neutral shell charge distribution and (ii) the 
same as (i) but with an added homogeneous background 
charge. Thus, when we will speak of a minimum energy 
it is assumed to be with respect to the shells' orientations 
unless stated otherwise. 
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Determining the orientation where the interaction en- 
ergy has a minimum can be numerically demanding since 
the distributions have a varying number of local extrema 
the number of which increases with I, and that also be- 
come shallower at larger intershell separations. This in 
general increases the probability of finding only the lo- 
cal minimum of the interaction instead of the global one. 
However, we have implemented a minimization procedure 
based on simulated annealinjpESI which reliably finds the 
global minimum in all of the cases we have considered, 
with an error in the interaction energy always below ap- 
proximately 5%. 



1. Neutral shells 

Considering the case of two neutral shells we have only 
one term in the interaction energy, Vu, which depends on 
several parameters: the shells' orientations, their sepa- 
ration, as well as the wave number I, the corresponding 
symmetry, and the strength of the variation <ii (equal for 
both shells) . Also variable are the bulk salt concentration 
Co and the radii of the shells R. 

To begin with, we show in Fig. [7] the absolute values 
of the minimum energy for the two shells at contact (p — 
2R); the interaction in the minimum is always attractive. 
The strength of the variation influences the scale of the 
energy in a predictable way, scaling with the square of ai 
[Eq. pi])]. Lowering kR increases the interaction energy, 
as already observed in Sec. |III[ and the energy decreases 
with increasing I. Also, at a fixed / the attraction in the 
minimum is stronger for a lower symmetry. 

The I dependence can be understood by considering 
the factors contributing to the interaction energy Vu in 
Eq. ( |31[ ), coming mainly from the function C(l,l, kR) 
and the modified spherical Bessel functions of the second 
kind k s (np), where the sum over s is limited by I (Ap- 
pendix All. The function k s takes on the biggest value 



when s = /, and thus the scaling of the energy is approx- 
imately set by the product of C(l,l, kR) and ki(Kp); the 
product which decreases as I increases. 

The symmetry dependence and to some extent the I 
dependence as well are likely linked to the average size 
of the local charge patches with the same sign. Lower 
symmetries at a fixed I have less minima, and the local 
patches of charge thus span a greater fraction of the shell 
surface. The interaction of two large patches brought 
close is less affected by neighboring patches some distance 
away, whereas the smaller patches of higher symmetries 
see an interaction averaged-out over several other patches 
in proximity. 

Figures [H] and [9] present the dependence of the interac- 
tion energy on the intershell separation for a number of 
different cases. The configuration of the two shells with 
the minimum energy does not appear to change with in- 
creasing separation; however, for every case there are sev- 
eral equal solutions due to the symmetry of the problem. 
We can observe again that lower I and lower symmetries 
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Figure 7. Absolute value of the minimum interaction en- 
ergy Vu at contact for two uncharged shells with symmetric 
variation (the interaction in the minimum is always attrac- 
tive). Circles show wave numbers with icosahedral symmetry, 
squares ones with octahedral, and diamonds those with tetra- 
hedral symmetry. Full symbols show the case where kR S3 10 
(co = 100 mM and R = 10 nm), and empty symbols the case 
with kR « 3 (co = 10 mM and R = 10 nm). Strength of the 
variation is taken to be er; = 1 eo/nm 2 . 



result in bigger attraction. 

From Fig. [8] it is obvious that the scaling of the interac- 
tion energy with the intershell separation at a fixed wave 
number is influenced only by the dimensionless param- 
eter kR, with the energy falling off more rapidly with 
larger kR. Changing the wave number also affects the 
scaling, as can be seen in Fig. |9j and the energy again 
falls off more rapidly with increased I. 

Since the depths of minima and maxima of the surface 
charge distributions coincide for odd wave numbers and 
differ for even wave numbers, something similar is also 
observed in the interaction free energy when comparing 
different shell orientations. For odd Is the interaction en- 
ergies in the minimum and maximum configuration mir- 
ror each other, even with changing separation, while for 
even Is the maximum of the interaction is usually some- 
what larger than the minimum. For all the other orien- 
tations of the shells the energy falls in between the limits 
delineated by the minimal and maximal interaction, all 
the while keeping a monotonic behavior with respect to 
the intershell separation, i.e. being either purely repul- 
sive or purely attractive. 
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Figure 8. Minimum of the interaction energy Va as a function 
of the intershell separation for symmetries with even wave 
numbers: icosahedral symmetry with I — 6 (circle) and / = 10 
(diamond) and octahedral symmetry with / = 6 (square) and 
I — 10 (triangle). Full lines show the case where kR ~ 10 and 
dashed lines a case with nR w 3. The inset shows the behavior 
of the absolute value of the energy for the two symmetries with 
I — 6, showing the scaling at larger separations. Strength of 
the variation is taken to be <r; = 1 eo/nm 2 . 



Figure 9. Minimum of the interaction energy Va as a func- 
tion of the intershell separation for symmetries with odd 
wave numbers: icosahedral symmetry with I = 15 (circle) 
and / = 21 (diamond) and octahedral symmetry with I — 9 
(square) and I — 13 (triangle). Full lines show the case where 
k,R m 10 and dashed lines a case with nR w 3. The inset 
shows the behavior of the absolute value of the energy for 
/ = 9 and I — 15, showing the scaling at larger separations. 
Strength of the variation is taken to be <t; = 1 eo/nm 2 . 



2. Charged shells 

On addition of a homogeneous background charge we 
can study the interaction of two charged shells in the 
presence of local charge density variation. The shells are 
assumed equal in size and the interaction energy is then 

VW(p, = V (p) + Wio(p,Ui) + Vu(p,u>i), (44) 

where I is again chosen from the allowed wave numbers 
for a given symmetry. In addition to the parameters 
considered in the case of two neutral shells we can now 
modify the total charge on the shells by changing <r , and 
we can assume that the ratio ai / ctq will play an important 



is expectedly bigger than the critical ratio in Eq. (40) 



role, similar to what we found out in Sec. Ill 



Taking a look first at the two shells at contact (p = 
2R), we show in Fig. 10 the contours of the minimum in- 



teraction energy as a function of kR and the ratio ai / <7q 
for the distributions already depicted in Fig. [5] Several 
things are notable: A sufficiently large charge variation 
with respect to <7o will cause a shift from repulsion to the 
attraction in the configuration of the two shells with the 
minimum energy. This critical ratio where the energy 
changes sign depends on kR, and becomes increasingly 
large when kR approaches 0. On the other hand, for large 
kR the critical ratio seems to approach a constant, which 



where the sign of the charge is still equal everywhere on 
the shells. Additionally, as already observed in previous 
cases, the sign of ai has an effect only on distributions 
with even wave numbers, and the value of cr simply sets 
the energy scale. However, the contours where the inter- 
action energy changes sign remain the same regardless of 
the value of <7o ■ 

The major influence on the critical ratio where the 
minimum of the interaction energy changes sign appears 
again to be the number of cxtrema in the distribution, 
where lower number of extrema needs lesser variation 
strength to achieve attraction - meaning lower I and 
lower symmetries. In distributions with the same sym- 
metry and similar / the minima of those with odd wave 
numbers are slightly stronger compared to distributions 
with even I. This can be understood by taking a look at 
the cross-term Wiq, where one can observe from Eq. (31 ) 
that a factor of (—1)' appears in the summation over the 
expansion coefficients. Thus, in the minimum configura- 
tions of distributions with odd I the cross-term is several 
orders of magnitude lower than the Vq and Vu contribu- 
tions, whereas the cross-term is comparable if the / is 



even (Fig. 11 ) 



From this Section we can conclude that even like- 
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Figure 10. Contour plots of the minimum of the interaction energy Vint in units of fcsT as a function of nR and the ratio 
(Ji/aa. The charge is kept fixed at erg = 0.4 eo/nm 2 , and we vary the variation ai. The shells are in contact, p = 2R. Dashed 
lines show the contours where the interaction changes sign, and dot-dashed lines show the critical variation ai [Eq. ( |40[ )] . Below 
the critical variation the charge on the shells is everywhere of equal sign, irrespective of the variation. The wave numbers and 
symmetries are the same as in Fig. [5] a) icosahedral with I = 15, b) octahedral with I = 9, c) icosahedral with I — 21, and d) 
dodecahedral with I — 6. 



charged particles can be found in configurations where 
they attract each other, if only the local variation of 
charge is large enough. It is the presence of patches of 
opposite sign that mediates this attraction, as can be 
seen from the regions excluded by dot-dashed lines in 
Fig. [TO} The attraction sets in for close separations of 
the two particles, as the long-range behavior is domi- 
nated by the repulsion of the homogeneous contribution, 
however small (Fig. 11). Any configuration other than 
the one with minimum energy can show either a purely 
repulsive or a mixed short-range attractive/long-range 
repulsive behavior, as already observed in Sec. |III| 



V. DISCUSSION AND CONCLUSIONS 

We have derived an analytical expression for the in- 
teraction of two permeable, arbitrarily charged spherical 
particles in the DH approximation. To study the inter- 
actions of virus-like particles or other shells with sym- 
metric variations in surface charge distributions we have 
introduced model charge distributions with icosahedral, 
octahedral, and tetrahedral symmetry, and characterized 
by a single wave number. 

We find that local variations in the charge distri- 
butions leading to local charge separation into oppo- 
sitely charged patches reduce the repulsion between like- 
charged shells, a feature that has also been observed 
in models of protein-p rotein interactions^ as well as in 
DNA-DNA interaction ;] 49 * 55 !. An important consequence 
of the reduced repulsion is that it can be furthermore 
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p/2R p/2R p/2R 



Figure 11. Total interaction energy Vint of two charged shells with icosahedral variation (full line) and the relative contributions 
(Vb - dotted line, Wio - dot-dashed line, Vu - dashed line) as a function of the intershell separation for a) I = 6 and ai/cro = 6, 
b) I = 15 and ui/ao = 6, and c) I = 15 and ai/ao — 10. The radii of the shells are R = 10 nm with k_R « 10, and the charge 
on the shells is fixed to ctq = 0.4 eo/nm 2 , making the Vb contribution equal for all three cases. 



turned into explicit attraction if the background homo- 
geneous charge is small enough. This onset of attraction 
between like-charged particles is due to the presence of lo- 
cal variation of charge, since the attracting patches have 
to have opposite sign. The attraction is short-ranged 
with respect to the separation between the particles, with 
the long-ranged interaction engendered by the homoge- 
neous charge background still being (vanishingly) repul- 
sive. Apart from the minimum interaction free energy 
configuration, any other configuration can exhibit either 
a purely repulsive or a mixed short-range attractive/long- 
range repulsive behavior. 

One of the drawbacks of the approach used here is the 
linearization of the non-linear Poisson-Boltzmann (PB) 
equation. The approximation is well-justified for radii 
relevant for viruses^ (R > 10 nm) and salt concentrations 
of Co ~ 100 mM. However, for smaller salt concentrations 
where kR ~ 1 the DH theory overestimate s the values of 
the potential and consequently the energjl2l2SI meaning 
that the values obtained for the interaction energy at 
shell contact are too big. The qualitative behavior, on 
the other hand, is expected to be much the sam^H. 

The non-homogeneous charge distributions introduced 
in this paper were all mapped on a single-shell capsid 
model. A further improvement would be a two-shell 
model, with each shell carrying a distribution corre- 
sponding to the charge on the inner and outer shells 7 . 
Nonetheless, the single-shell model is able to capture the 
main aspects of the electrostatic interactions in the sys- 



tenP 1 . 

If the charge distributions can be characterized by 
one or only a few different wave numbers (like the ones 
considered in this paper), this greatly simplifies the ex- 
pressions for the interaction potential, limiting many of 
the sums involved. On the other hand, distributions 
where patches extend over well-defined surface regions 
can be difficult to be parametrized within the formalism. 
The most extreme cases would be distributions of point 
charges on a shell or a partially formed capsid as con- 
sidered in Ref. 3 - i.e. distributions consisting of delta 
functions or Heaviside step functions in the solid angle, 
which are also notoriously difficult to accurately repre- 
sent with a finite number of terms in an expansion. 

The symmetric charge distributions presented can 
serve as improved models in a broad variety of virus 
systems. One such example could be the adsorption of 
flexible polyelectrolytes on viral capsids which have do- 
decahedral charge distribution, where a combination of 
the lowest dodecahedral symmetry (1 — 6 in Fig. [5| with 
the addition of a homogeneous background nicely fits the 
model of discrete charges as used in Ref. [5(3 These mod- 
els also open up the possibilities of studying the inter- 
action of viral capsids with multivalent ions within the 
dressed counterions theory, or even multivalent ion me- 
diated interaction between two such capsids^HSSl. 

The derived expressions for the interaction potential 
are certainly too heavy to be used in computer simula- 
tions. And even though molecular dynamics simulations 
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show that switching from continuous to discrete charge 
patterns in model spherical proteins leads to significant 
differences between the resulting screened Coulomb in- 
teractions^, the continuous distributions presented are 
more tractable analytically, especially since the symme- 
tries can be described by using one or at most few differ- 
ent wave numbers. 
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Appendix A: Details of the derivation 



1. Properties of Wigner 3-j symbols 



The combination of the two Wigner 3-j symbols in the 



function Lf^ [Eq. ( 22 )] gives the rise to following proper- 



tie d 37 * 38 !: the function is non-zero only when m = q and 
I + p + s = even. The product of the two Wigner 3-j 
symbols also makes the function invariant under all 
odd and even permutations of the three wave numbers 
I, s, and p. The sum over s is limited with the triangle 
inequality, \l — s\ < p < I + s. This further restricts s to 
even numbers when / = p; on the other hand, when one 
of them is zero (e.g. p = 0), the sum contains only one 
term, s = I. 



2. Obtaining coefficients for the electrostatic potential of 
the shells 



Writing out the continuity equation [Eq. (15l] on the 



surface of the second shell in terms of the functions ip lm 
we obtain 

^2a(lm,2)ii(KR)Yi m {Q) = 

l,m 

= J2Hlm,2)k l { K R)Y lm (n) + 

l,m 

+ ]T b(im, i) ]T £f£(i); P M)y M (o). (Ai) 



By multiplying both sides with Y^ m , and using the or- 
thogonality properties of spherical harmonics^, we have 
after integrating over Q 



a(lm, 2) - J2 Hpq, 1)^7(1) U{kR) = b(lm, 2)ki(nR). 
p,q 

With the same procedure the boundary condition for the 
discontinuity of the electric field [Eq. ( 16 1] on the surface 
of the second shell on the other hand yields 



2) -X>(pg,i)/:£(i; 



i'i(nR) 



-bilm^k'^KR) 



<j(7m,2) 

KE£q 



(A3) 



The factors in the square braces of Eqs. (A2) and (A3) 
are the same, giving 

.kiinR) ., „. , .... „, ailing) 
b(lm,2) /; ( i' l (KR)-b(lm,2)k' l (KR) - 



ii(nR) 



K£Eq 



(A4) 

from which we obtain the coefficients b(lm, 2) of the ex- 
terior solution [Eq. (23)]. Same approach is used to ob- 
tain b(lm,\). By inserting these coefficients back into 
the boundary condition equations, we then get the coef- 
ficients of the interior solution, Eq. (24). 



3. Deriving the symmetrized terms of interaction energy 



The function Wi p we have introduced in Eq. (31 1 written out in full is 



2nesr i 

ra,q \ 

Cq{1, KR)ip(nR) 
ki(nR) 



k p (KR) 



a(pq,l)a*(lm,2)£ l ™(l) + a*(lm,l)a(pq,2)£ l ™(2) 



a(lm, l)a*(pq, 2)£f« (1) + a*(pq, l)a(lm, 2)£f« (2) 



(A5) 
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We can separate the four terms of the sum into two parts, Wi p = W l 



(i) 
i P 



w[p\ where 



< 3 = 



2neeo 



E 



R 2 ^ [ C (j>, kR)h(kR) o\ri™fT\ , C (l,nR)i p (KR) 
k j^-ftj cr{pq,l)a (lm,2)£ pq (l) + 



a>g,l>(Zm,2)£^(2) 



(A6) 



and similarly for . The Wigner 3-j symbols appearing in the functions £^ are non-zero only if q = m. Since 

the functions are invariant under odd and even permutations we get that (1) and £f™(2) differ only in the factor 
(— l) s appearing in the sum in the latter case. Thus the sum over q disappears and using the fact that l+p + s = even, 
we get 



W, 



R 2 



J2(-iy+ m (2s + l)v/(2/ + l)(2p+l)fc s M 



2nee 
C (p, KR)ii(nR) 





s 













! 




m y 



/ -t \ */, 0A . Ca(l,nR)i p (KR) 
<T(pm,l)a (lm, 2) H — r — <r (pm, l)a(lm, 2) 



k p {nR) ' ki(nR) 

Additionally we can use the following relations: 

Co(p, nR)ii(nR) Cq{1, KR)ip{nR) 

kp(nR) ki(nR) 

Co(p, KR)ii(nR) Co(l , KR)ip(KR) 

kp(nR) ki(nR) 

where the function C(l,p,x) is defined in Eq. ( |32[ ), and after some manipulation we obtain 

R 2 



2C(l,p, kR), 
0, 



(A7) 

(A8) 
(A9) 



KEEq 



C(l,p, kR) ^(-l) p+m Re[cr(pm, l)a*(lm, 2)] 



x (2s + l)V(2J + l)(2p+l)fc a («p) ( o ) ( 



/ s p 
m m 



(A10) 



(2) 

In an analogous manner we can derive a similar expression for , and putting both results together we retrieve 



Eq. (31) 



Appendix B: Limiting cases of interaction energy 

In this Appendix we shall derive some limiting cases of 
the general expression for the interaction energy to show 
that we obtain the results found in the literature. A 
useful partial result is the expression for the interaction 
of two homogeneously charged shells, Vq. 

The distribution is spherically symmetric and invariant 
under all rotations, thus giving for both shells a(00,i) = 
\[A~koq (this normalization of the coefficient is due to the 
requirement that AttR 2 ^ — Q gives the total charge on 
a shell). Taking this into account, the interaction of two 
such shells is found to be 



TA . , 4nR 2 a 2 . , 2 . e" K " 
Vo(p) = -smh'UR) x . 

K££0 Hp 



(Bl) 



as 



lim C(l,p, x) 

i<0 



J+p_ 



2 -(/+p+l) a .2 



r(z + 3/2)r(p + 3/2)' 



(B2) 



where T(z) is the Gamma function^. The bigger the 
wave number of the charge variation, the smaller role it 
plays in this limit, and in the lowest order of the expan- 
sion in x = kR we need focus only on the lowest wave 



numbers I = p = 0. Thus we obtain using Eq. (Bl ) 



lim V(p,LJi 

KR->0 



Q2 

lim Vo(p) = 



P 



(B3) 



the screened Coulomb (DH) interaction of two point 
charges Q. 



1. Coulomb limit 



2. Limit of large screening 



Firstly, we take a look at the Coulomb limit where 
kR — > 0. In this limit, the function C(l,p,x) found in 
partial terms of the interaction energy Wi p [Eq. pi])] goes 



Another limit is the regime of large screening, kR. — > 
oo. We need to consider the asymptotic expansions of 
the function C(l,p,x) and the modified spherical Bessel 
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functions appearing in the expression for the interaction 
energy: 



(_l)' + (-lV 



2tt 



2tt 



k s (Kp) 



tt e 



2kH + 2R 



4) 



(B5) 



where we have written the distance between the shell 
centers p in terms of their distance of closest contact 
H : p — H + 2R. From Eq. (B5) it is obvious that in 



this limit the dominant term in the radial dependence is 
independent of the wave number, but the anisotropy of 
the two distributions persists-^. This makes it hard to 
obtain any general results as the interaction in general 
depends on the orientations of the two shells when they 
are brought together. But if we consider two shells at 
small separations, p <C i?, we can check the interaction of 
two homogeneously charged shells [Eq. (Bl|] and obtain 



V (H) = 



8/v 2 ££o 



(B6) 



which has the same functional dependence as given 
by Verwey and Overbeek^, but differs in the prefactor 
due to a different approximation. 



3. Dipole-dipole interaction 

Another limit we can check is the expression for the in- 
teraction of two (equal) dipole distributions. We concern 
ourselves with an axially symmetric charge distribution 
with I = 1, cr(lm, i) — aid m Q. The connecting line be- 
tween the two dipoles is chosen to be z, and due to the 
symmetry of the problem the only relevant Euler angles 
are the azimuthal angles of the two dipoles, fa and fa. 
The dipole moment of the shells can be calculated from 
their surface charge distributions^, 



ra(n)6(r-R) d 3 



(B7) 



The distribution in the reference frame is axially sym- 
metric, and we get a dipole moment in the z direction 
only with the magnitude 



fi = aiR 3 



(B8) 



Thus, we can write for the general orientations of the 
two dipoles ^i/x 2 = p, 2 (cos fa cos fa + sin fa sin fa ) and 
(/xiz)(/x 2 z) = p? cos fa cos fa • After summing over m 
and s the dipole-dipole interaction energy is 

p2 2 

1 C(1,1,kR) 



V n = 



Sire- Kp 



kesq ' ' \ 2(np) 3 

x (2 + 2np + (up) 2 ) cos fa cos # 

- (1 + Kp) sin fa sin fa 



(B9) 



Simplifying the above equation using the expressions for 
the dipole moment we can write in the Coulomb limit of 
two point dipoles 



lim Vn/k B T = l B - 



- Kp 



kR^O 



(1 + K,p)ni(J,2 



- (3 + 3k P + M^C/tix^^z^BlO) 

obtaining the same expression for the DH interaction of 
two point dipoles as given in the literature (e.g. Ref . IS0| . 
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